Journal of Quantitative Spectroscopy & Radiative Transfer 113 (2012) 2351-2355 



ELSEVIER 


Contents lists available at SciVerse ScienceDirect 

Journal of Quantitative Spectroscopy & 
Radiative Transfer 

journal homepage: www.elsevier.com/locate/jqsrt 


a! 

ournal of 
uantitative 
pectroscopy & 
adiativc 
ransfer 


Scattering and absorption properties of polydisperse 
wavelength-sized particles covered with much smaller grains 

Janna M. Dlugach 3 , Michael I. Mishchenko b '*, Daniel W. Mackowski c 

a Main Astronomical Observatory of the National Academy of Sciences of Ukraine, 27 Zabolotny Street, 03680 Kyiv, Ukraine 
b NASA Goddard Institute for Space Studies, 2880 Broadway, New York, NY 1 0025, USA 
c Department of Mechanical Engineering, Auburn University, AL 36849, USA 


ARTICLE INFO 


ABSTRACT 


Available online 17 March 2012 

Keywords: 

Maxwell equations 
Electromagnetic scattering 
Polarization 

Scattering and absorption characteristics 
Small-scale surface roughness 
Wavelength-sized particles 


Using the results of direct, numerically exact computer solutions of the Maxwell 
equations, we analyze scattering and absorption characteristics of polydisperse com- 
pound particles in the form of wavelength-sized spheres covered with a large number 
of much smaller spherical grains. The results pertain to the complex refractive indices 
1.55 + i0.0003, 1.55 + i0.3, and 3 + iO.l. We show that the optical effects of “dusting” 
wavelength-sized hosts by microscopic grains can vary depending on the number and 
size of the grains as well as on the complex refractive index. Our computations also 
demonstrate the high efficiency of the new superposition T-matrix code developed for 
use on distributed memory computer clusters. 

Published by Elsevier Ltd. 


1. Introduction 

The growing efficiency of computers and numerically 
exact codes based on direct solutions of the Maxwell 
equations allows one to study scattering and absorption 
properties of small particles with increasingly complex 
morphologies [1-11]. In recent papers [12,13], we intro- 
duced a new particle morphology in the form of a wave- 
length-sized spherical host covered with a large number of 
much smaller spherical grains. This morphology can be 
used to study the potential optical effects of microscopic 
surface irregularities of wavelength-sized particles 
(cf. [14-24]) as well as to model the actual “dusting” of 
wavelength-sized particles by microscopic grains [13], 

The analysis in Ref. [12] was limited to one refractive 
index typical of mineral aerosols (m = 1.55 + i0.0003, 
where i =(— 1 ) 1/2 ) and to fewer than 50 small grains 
covering a large host. This paper is a natural extension 
of Ref. [12] in two respects. First, we analyze the effect of 
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a five-fold increase in the number of small grains, where- 
upon the grains cover the surface of the large hosts almost 
entirely. Second, we analyze the consequences of increas- 
ing absorption (m = 1.55 + i0.3) as well as increasing both 
the real and the imaginary parts of the refractive index 
(m = 3 + iO.l). A related objective is to further assess 
the efficiency and versatility of the new public-domain 
Fortran-90 code based on the superposition T-matrix 
method (STMM) [25] and developed for use on distributed 
memory computer clusters [26,27], 

2. Numerically exact computer modeling 

Following Ref. [12], we consider a scattering target 
obtained by placing a number N g of small non-overlap- 
ping spherical grains of radius r randomly on the surface 
of a relatively large spherical host of radius R using a 
random number generator (see the inset in Fig. 1(a)). The 
refractive indices of the host and the grains are the same. 
It is assumed that the target is illuminated by a parallel 
quasi-monochromatic beam of light and that the observa- 
tion point is located in the far-field zone. It is also 
assumed that all scattering and absorption characteristics 
are averaged over the uniform orientation distribution of 
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the resulting host-grains configuration with respect to 
the laboratory coordinate system. The transformation of 
the Stokes parameters upon far-field electromagnetic 
scattering is then written in terms of the normalized 
Stokes scattering matrix [28]: 


- jsca - 


arid) b i(0) 0 0 


jinc 

Q sca 

U sa 

oc 

fq(0) a 2 (0) 0 0 

0 0 a 3 (0) b 2 (0) 


Q inc 

U inc 

V sca 


0 0 -b 2 (0) a 4 (0) 


v inc 


( 1 ) 


where 0e[O°,18O°] is the angle between the incidence and 
scattering directions, while both sets of the Stokes para- 
meters are defined with respect to the scattering plane. 
Strictly speaking, the off-block diagonal elements of the 
scattering matrix vanish only if the compound particle 
has a plane of symmetry. We have found, however, that in 
the cases considered, the elements denoted by zeros are 
negligibly small (in the absolute-value sense) in compar- 
ison with the other elements at the same scattering 
angles. The (1,1) element a^O), called the phase function, 
is normalized according to the integral condition 

l- f a 1 (0)sin0d0= 1. (2) 

2 Jo 


This model of compound aerosols allows us to 
compute all scattering and absorption characteristics 
using the highly efficient and numerically exact STMM. 
Specifically, we used the recently developed parallelized 
computer program posted on-line [27] to perform exten- 
sive computations for different values of JV g and m. Half of 
our results pertain to the refractive index 1.55 + i0.0003 
typical of mineral aerosols in the visible part of the 
spectrum. In order to analyze the effect of increasing N g , 
we consider the values N g = 0, 49, 149, and 249. To study 
the effect of increasing absorption, computations were 
also performed for m = 1.55 + i0.3. Finally, we consider 
the case of hematite particles with m=3 + i0.1 (cf. 
Refs. [21,23]). In the majority of simulations, the ratio of 
the host radius to that of the small grains was kept constant 
at R/r= 10. However, to analyze the effect of increasing grain 
size relative to that of the host, we also considered compound 
particles with m= 1.55 + i0.0003 and R/r= 5. The wavelength 
was kept constant at 0.62832 pm. 

To eliminate resonances typical of monodisperse par- 
ticles [16,28], all scattering and absorption characteristics 
were averaged over the standard power-law distribution 
of host radii [28]: 


n(R) = 


constant x R 3 , R min < R < R max 
0 otherwise 


(3) 


with an effective radius of R e ff=l pm and an effective 
variance of v eff =0.05. Of course, averaging over sizes also 
renders a better representation of polydisperse particle 
ensembles encountered in various applications. The aver- 
aging was based on the Gauss-Legendre quadrature 
formula with 100 division points applied to the corre- 
sponding range of host radii 0.661 pm < R < 1.439 pm. 
The corresponding effective radius of the small grains 
(equal to 1/5 or 1/10 times R e ff) will be denoted by r e ff. 


3. Numerical results and discussion 

Table 1 lists ensemble-averaged integral radiometric 
characteristics, viz., the extinction, C ext , scattering, C sca , and 
absorption, C abs , cross sections as well as the single-scatter- 
ing albedo vj and the asymmetry parameter g. The data 
presented in the first four rows of the table show how 
these characteristics change with N g in the case of 
m = 1.55 + i0.0003. It is seen that all three cross sections 
increase, especially Qxt and C sca , while the single-scattering 
albedo and the asymmetry parameter hardly change. The 
significant increase of the extinction and scattering cross 
sections can be attributed primarily to the growing geome- 
trical cross sections of the compound particles. Increasing 
absorption in both the hosts and the small grains 
(m=1.55 + i0.3, JV s =0 and 49) results in a non-negligible 
increase of all three cross sections and a decrease of the 
single-scattering albedo with increasing N g (compare the 
sixth and seventh rows). Comparison of the first, second, 
and fifth rows reveals a much stronger effect of dusting the 
large hosts with R/r= 5 grains. All three cross sections 
increase quite substantially, whereas the asymmetry para- 
meter is significantly reduced. Finally, the last two rows of 
Table 1 reveal rather modest changes in the integral radio- 
metric characteristics upon covering R e ff= 1 hematite hosts 
with N g =49 small (r/R=0.1) hematite grains. 

Fig. 1 parallels Table 1 and depicts non-zero elements 
of the Stokes scattering matrix for both “clean” and 
“contaminated” hosts. Fig. 1(a) shows that the effect of 
increasing N g on the phase function is three-fold: (i) the 
forward-scattering diffraction peak becomes somewhat 
stronger; (ii) there is a modest increase at side-scattering 
angles; and (iii) the backscattering peak is suppressed 
substantially, the differences between the N g = 0 and 
N g =249 phase functions exceeding a factor of 4. By the 
time N g reaches the value 249, the small grains form a 
quasi-uniform coating around the hosts, thereby making 
the compound particles quasi-spherical with R e ff ~ 1-2 pm. 
However, the resulting increase in the particles’ geome- 
trical cross sections can explain only trend (i), as the gray 
curve in the upper left panel of Fig. 1(a) vividly demonstrates. 
Trends (ii) and (iii) are, in fact, reminiscent of the typical 
effects of increasing asphericity of such morphologically 
simple particles as spheroids [28-30]. The phase-function 
results computed for the case of R/r= 5 and N g =49 are even 
more indicative of these effects (see the dark blue cun/e in 
the upper left panel of Fig. 1(a)), which is not surprising given 
the substantially larger size of the grains leading to visibly 
more pronounced asphericity of the resulting compounded 
particles. The linear-polarization curves in Fig. 1(a) reveal 
another typical effect of increasing asphericity resulting in 
progressively neutral polarization [31], The - b ] (0)la ] (J)) 
curve plotted for R/r=5 and N g = 49 even appears to show 
the gradual development of a positive-polarization maximum 
at side- and near-backscattering angles typical of spheroids 
and natural mineral aerosols [31-33]. 

The upper right-hand panel of Fig. 1(a) depicts the 
computed ratio a 2 (0)/a 1 (0), which is identically equal to 
100% for “clean” spherical hosts but deviates from 100% 
for contaminated particles, thereby signifying the non- 
sphericity of the latter and causing non-zero linear 
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Scattering angle (deg) Scattering angle (deg) 

Fig. 1 . (a) Elements of the normalized Stokes scattering matrix for "clean" polydisperse hosts and polydisperse hosts "dusted" with different numbers of 
small grains. The refractive index is 1.55 + i0.0003. The inset illustrates the morphology of the compound particles for N g =49. (b) As in (a), but for 
refractive indices 1.55 + i0.3 and 3 + iO.l. The effective radii of the hosts and the small grains are fixed at 1 pm and 0.1 pm, respectively. 


depolarization ratios [34]. The deviation of the ratio a 2 (0)/ 
ai(0) from 100% is especially large in the case of R/r= 5 
and JV g =49, which can be explained qualitatively as being 
the consequence of the greater overall asphericity of the 
resulting compound particles and stronger “electromag- 
netic interactions” between the grains and the host. 


The most obvious manifestation of the nonsphericity 
of the compound particles in the a 4 (0)/ar(0) curves is their 
deviation from -100% at 0=180°, causing non-zero 
circular depolarization ratios [34], Again, this deviation 
is especially large for compound particles with R/r= 5 and 
JV g =49. On the other hand, the numerical data depicted in 
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Table 1 

Ensemble-averaged extinction, scattering, and absorption cross sections per compound aerosol particle, single-scattering albedo ro, and asymmetry 
parameter g as functions of the number N g of small grains covering the surface of a large host. 


m 

feff/Keff 

N g 

C ext (pm 2 ) 

C sca (pm 2 ) 

Cabs (pm 2 ) 

zu 

g 

1.55 + i0.0003 

_ 

0 

6.314 

6.265 

0.049 

0.9923 

0.668 

1.55 + i0.0003 

0.1 

49 

6.607 

6.558 

0.049 

0.9926 

0.661 

1.55 + i0.0003 

0.1 

149 

7.358 

7.304 

0.054 

0.9927 

0.661 

1.55 + i0.0003 

0.1 

249 

8.214 

8.156 

0.058 

0.9931 

0.669 

1.55 + i0.0003 

0.2 

49 

10.670 

10.606 

0.064 

0.9940 

0.609 

1.55 + i0.3 

- 

0 

6.464 

3.183 

3.281 

0.492 

0.912 

1.55 + i0.3 

0.1 

49 

6.845 

3.265 

3.580 

0.477 

0.920 

3 + iO.l 

- 

0 

7.355 

4.312 

3.043 

0.586 

0.776 

3+iO.l 

0.1 

49 

7.448 

4.408 

3.040 

0.592 

0.736 


the bottom right-hand panel of Fig. 1(a) demonstrate 
rather insignificant changes of the ratio b 2 (d)la^d) with 
increasing JV g and/or increasing grain size. 

In order to analyze the effect of increasing absorption, 
in Fig. 1(b) we present the numerical data obtained for 
R e ff=l pm, R/r= 10, JV g =0 and 49, and m = 1.55 + i0.3. 
Comparison with the analogous cases in Fig. 1(a) reveals a 
similar magnitude of changes caused by the presence of 
N g =49 grains on the surfaces of much larger hosts. 

Finally, the combined effects of increasing both the 
real and the imaginary parts of the refractive index are 
revealed by the results of computations for “clean” 
and “contaminated” hematite particles with R e ff=l pm, 
R/r= 10, N g =0 and 49, and m = 3 + i0.1 (Fig. 2(b)). In 
contrast to the integral radiometric characteristics dis- 
cussed above, dusting each large host with 49 small grains 
causes a noticeable change in all elements of the scatter- 
ing matrix. The deviation of the ratio a 2 (d)la^(d) from 
100% at side-scattering angles is especially pronounced. 

4. Concluding remarks 

This paper in combination with Refs. [12,13] demon- 
strates that large host particles covered with a large 
number of microscopic grains represent an interesting 
morphological class of wavelength-sized scatterers which 
can be used to model real “dusting” of larger aerosols by 
much smaller ones or, possibly, the optical effects of 
microscopic surface irregularities. In the case of a modest 
number of grains and a small real part of the refractive 
index, the optical consequences of such “contamination” 
are relatively weak. They become more pronounced with 
increasing real part of the refractive index, with increasing 
grain size relative to the host’s size, and/or in cases when 
the smaller grains cover a large fraction of the host’s 
surface. In some respects the optical effects of increasing 
the number of grains or the r/R ratio in the case of the 
refractive index 1.55 + i0.0003 appear to be qualitatively 
similar to those of increasing the aspect ratio of spheroids, 
although a quantitative analysis of this similarity is 
problematic given the obvious geometrical incompatibility 
of these two types of particle morphology. By far the greatest 
are the effects of surface contamination of weakly absorbing 
mineral hosts by strongly absorbing soot grains [13], 

A natural byproduct of our study was a further 
assessment of the utility of the new Fortran-90 STMM 
code [26,27] when run on distributed memory computer 


clusters. Based on numerous timing tests, we conclude 
that the requisite computer time is inversely proportional 
to the number of nodes and to the number of cores in 
each node. Overall, our extensive experience with the old 
and new versions of the STMM code demonstrates that 
the latter is both highly accurate and very efficient. 
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